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Abstract 
> 

A theoretical description of the phenomenon of modulation of near- 
*"** ■ wall turbulence by large scale structures is investigated. The descrip- 

' . ' tion given is simple in that the effect of large-scale structures is limited 

m 

o 

£\j . scale fluctuations of the skin friction. The most natural and compact 

form of expressing this mechanism is given by the usual Reynolds- 



>< 



c3 



number-independent representation of the total skin friction and ve- 
locity, scaled in wall variables, where the mean quantities are replaced 
jj| , by large-scale low-pass-filtered fluctuating components. The theory is 



rewritten in terms of fluctuations via a universal mean velocity and 
random mean square fluctuation velocity profiles of the small-scales 
and then linearised assuming that the large-scale fluctuations are small 
as compared to the mean components. This allows us to express the 
superposition and modulation coefficients of the empirical predictive 
models of the skin friction and streamwise fluctuating velocity given 
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respectively by Marusic et al. (13 i/l Eur. Turb. Conf., 2011) and 
Mathis et al. (J. Fluid Mech. 2011, vol. 681, pp. 537-566). It is 
found that the theoretical quantities agree well with experimentally 
determined coefficients. 



1 Introduction 



Direct numerical simulations of turbulent skin friction drag reduction are 
necessarily done at moderate Reynolds numbers (Re). Until recently, the 
hope that theories justified by comparisons with such calculations can be 
extrapolated to much larger Re typical for flight regime, was based on the 
idea of the universality of near- wall turbulence, where the skin friction is pri- 
marily generated. This hope suffered a serious blow with the discovery of the 
emergence, at sufficiently large Re, of a secondary peak in the pre-multiplied 
energy spectra map of t he streamwise velocity c ompon ent, namely the "outer- 
peak" (see for exa mple iHutchins fc Marusic I. l2007al). Empirical fits based 
on the latest data ( jAlfredsson. Orlii. fc Segalini 1120111 1 suggest that at high 
Re more turbulence kinetic energy is contained in large-scale structures than 
in the near-wall region where the inner-peak is located. In the past decade 
or so, the large-scale structures and t heir interaction with near-wall turbu- 
lence have attracted much attention (jAdrian. Meinhart. fc T omkins. 2000; 



del Alamo fc Jimenez ll2003uGanapathisubramani. Longmire. fc Marusic 
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2003; 



2006; 



Hambleton. Hutchins. fc Marusic I J2006I : iMcKeon fc Morrison 1 120071 : 



Morrison 1120071 . and even this list is far from complete). In particular, it has 
been shown that the large-scale structures associated with the log-layer sig- 
nific antly affect the near- wall tur bulence. lAbe. Kawamura. fc Choi I (120041 ) 
and IHutchins fc Marusic I (j2007al ) have observed that these large-scale mo- 
tions impose a substantial "footprint" all the way down to the wall, and this 
has been confirmed in recent direct numerical simulations (DNS) of evolving 
boundary layers (Orlii & Schlatter 2011). This footprint is seen as a long- 
wavelength co mponent superimposed onto the small-scale near-wall fluctua- 



tions. Further. IHutchins fc Marusid (1200761 ) and lMathis. Hutchins. &: Marusic 



(120091 ) demonstrated that the large-scale motions are not merely superim- 
posed near the wall, but also amplitude modulate the small-scale motions, 
and developed a tool to quantify the degree of modulation. Based on the 



above observations, Marusic and colleagues developed a simple quantitative 
model of the effect of large-scale structures onto the near-wall turbulence 
(Marusic et al. 2010; Mathis et al. 2011), and from here on we refer to these 
papers as MHM. In this model the effect of the large-scale structures is ex- 
pressed as a combination of an amplitude modulation and superposition of 
the large-scales onto the near-wall small-scales. The amplitude modulation 
factor and the superimposed additive term were approximated as linear func- 
tions of the low-pass-filtered velocity in the logarithmic layer. The model was 
shown to reconstruct a statistical representation of the near-wall fluctuating 
streamwise velocity component, as well as the fluctuating skin-friction, with 
good accuracy over a large range of Re. 

In the present paper we demonstrate that the modulation effect as ob- 
served in an experiment can be derived from a simple hypothesis. 



2 Skin friction modulation 

The classical view on the universality of near- wall turbulence consists in the 
statement that all the flow variables, if expressed in wall units, are indepen- 
dent of the Reynolds number. For skin friction r this means that 

- * (tv% xu T yu T \ 

T = TT , , , (1) 
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where r*(t + , x + , y + ) is a universal function of its arguments in the sense that 
all its statistics are independent of Re, f is the mean skin friction, u T = y/f/p, 
v is the kinematic viscosity, p is the density, and t, x, and y are time and 
coordinates in the plane of the wall. Note that all the quantities here are the 
total values: in this paper fluctuations will always be explicitly marked with 
prime or tilde. Note that the time average of r* equals one: r* = 1. 

The experimentally observed effect of large-scale structures on the near 
wall turbulence is in contradiction with (TTJ). We propose therefore to replace 
flTJ) with a formula recognising the dependence of the skin friction on the 
large scale effects. Our proposal might be reminiscent of the transition from 
the famous Kolmogorov-41 theory to Kolmogorov-62 theory in the isotropic 
turbulence, but we will not further comment on this similarity. 



The proposed hypothesis consists in the statement that the effect of large 
scale structures on the near-wall turbulence is limited to replacing f (and, 
correspondingly, u T ) in flTJ with slowly- varying Tr,{t) (and u TL (t)): 

r=rL(t)r ^^M,^w,^wy (2) 
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The slowly- varying component of the wall skin friction can be extracted from 
the full skin friction by a suitable filtering procedure. A quasi-steady effect 
described by (|2D can naturally be called a modulation of near- wall turbulence 
by large-scale structures. We will now try to verify the hypothesis (|2J) by 
comparing it with the available experimental data. 

Note that in Marusic et al. (2011) the modulation was understood some- 
what differently, and it was described in terms of fluctuations. To compare, 
we rewrite (J2J) in terms of fluctuations, with r* = 1 + r'* (for brevity we will 
omit arguments of r*). Assuming that tl and r* are not correlated, we also 
have r = T[. Then tl = r + r' L and r —f + r'. Substituting these formulae 
into (]2D arid rearranging gives 

/+ = ^= (l -4) r" + 4 (3) 



Compare this with equation (1) in Marusic et al. (2011), which has the 
form 



< 



(l + /Vj t ) / + au'+ OL , (4) 

where u' OL is the fluctuating large-scale signal from the log region, and we 
added primes to fluctuation quantities, to bring it closer to our notation 
(subscript "p" describing a predicted quantity). 

One can see that if our explanation of the modulation is true, then 



T, 



au . 



./+ 



' OL 
l 

and 

j = PW + OL , 



that is a and (3 should be equal. In Marusic et al. (2011) these quantities 
were found independently from experiment using very different procedures, 



and it was obtained that a = 0.0898 and (3 = 0.0867. The small difference 
between these values gives the first confirmation of our interpretation of the 
modulation effect. 



3 Velocity modulation 



The classical view on the universality of velocity distribution gives: 
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The idea of modulation can now be expressed again as a statement that 
the effect of large scale structures on the near-wall velocity distribution is 
limited to introducing into the above formula a dependence of f or, more 
conveniently in this case, the corresponding dependence of u T on time, with 
the related assumption that the variation of these quantities is much slower 
than the variation of u* : 



u\ * ( ^l^) xu TL (t) yu TL (t) zu TL (t) \ . . 

U — U TL [t)U I , , , \p) 
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Accordingly, this means that the large scale structure effect amounts to 
the amplitude modulation of the universal total (not fluctuating) velocity 
u* via the factor u T = u TL (t), a frequency modulation due to what in fact 
is the time-dependence of the wall units in the expression tu^. (t)/v, and a 
scale modulation due to the time- dependence of the wall units in xu TL (t)/u, 
yu TL (t)/u, and zu TL ii)lv. In treating the skin friction in Section |2l as far as 
the random mean square of the skin friction fluctuation was concerned, the 
frequency and scale modulation could be ignored because the skin friction 
random mean square is actually independent of time and spatial coordinates. 
In the present case, this is true only for the frequency and scale modulation 
in the wall-parallel directions. When, as a result of the variation of u TL the 
wall- normal non-dimensional distance zu TL (t)/u varies, so does the random 
mean square. For this reason in what follows we will frequently omit the 
tv%..(t)/v, xu TL {t)/v, and yu TL (t)/u arguments of u*, but we will often keep 
zu rL {t)/v. 



It is also convenient to introduce new variables reflecting the frequency 
and scale modulation, namely 

T+ = *<(*) z+ = zu TL (t) 

V V 

Note the difference between T + and Z + and the variables in wall units, 
defined similarly but with the time-averaged value of the friction velocity 
u T = u TL (t) : 



tut i zu 7 
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Our goal is to rewrite (JSJ) in fluctuations and to compare it with the 
formula given in MHM. 

Let us introduce the universal mean profile as 

U*{z + )= lim^- r T u*(t + ,z + )dt + (6) 

and the universal velocity fluctuation function 

u*(t+, z+) = u*(t + , z+) - U*(z + ). (7) 

Note that both the average and the fluctuation are introduced here with z + 
held constant as time varies. Where the velocity at a fixed distance z to the 
wall will be considered, the argument of u* will be Z + , which is varying in 
time, and not z + . This is why a tilde and not a prime is used to mark a 
fluctuation in (J2J). 

Note now the due to the separation of time scales of near-wall turbulence 
and large scale structures 



u*(T+,Z+) = U*(Z+). 

The assumption of the separation of time scales can be formalised and then 
the above statement can be proved formally, but in this paper we will fo- 
cus at the level of relating to the physical understanding. Introducing the 
fluctuation of the friction velocity u' = u TL (t) — u T , one can write 

u = (u T + u' TL (t)) (U* (Z + ) + u* (T+, Z+)) 
6 



Due to the separation of scales 



u' TL (t)u*(T+,Z- 
Hence the velocity average is 



u = u T U* (Z+) + u' TL U* (Z+). 
and its fluctuation, expressed in wall units, is 

u' + = V —^ = u* (T + , Z+) (l + ^£T\ + U* (Z + ) ^ + U* (Z+) - *- 

(8) 
The corresponding formula in MHM has the form 



/ + * 

u p — tt MHM 



[l + u u'i L }+a u u , i L , (9) 



in which we include primes where appropriate, and have added the subscript 
u to a and (3 to distinguish them from a and (3 of Section [2j 

Both these formulae represent the velocity fluctuation as a fast varying 
term multiplied by a slowly varying (depending only on large scale structure) 
amplitude modulation term plus a slowly varying superposition term. We 
might now try to associate u'(t)/u T with u'q L similar to Section [3 However, 
there are two difficulties here. First, in MHM (3 U was chosen in such a way 
that «y HM is not amplitude-modulated by the large scale structures. More 
precisely, (3 U is chosen in such a way as to minimize the effect of large scale 
structures on the random mean square of the fast- varying term. This can be 
reproduced in our formula by introducing the mean square of u* as 



u* rms (z + ) = ^ 2 (T+,zH 



and assuming that 



M MHM _ U rms( Z )~~2 ( 7+\ (• / 

Note that the right hand-side of fflOl) is frequency modulated via T + , while 
MHM did not mention frequency modulation. The same is true for scale 
modulation via xu TL {i)/v and yu TL (t)/v. This would make a difference where, 



for example, the spectra are concerned. However, for random mean square 
this difference does not matter. Then (JS} takes the form 

H- _ ..• Krns(Z + ) r , K L (t)\ rn u' TL (t) 



u p — M MHM7 



U rms\ Z ) 



(i + <M) + u* ( Z +) ^rM + u* t z +) _ 1 (ii) 



The second difficulty is that if we now try to relate the modulation factor 
and the superposition term to u'q L , at least one of them will have to depend 
on u'q L nonlinearly, while in (J5J) these relations are both linear. Noting that 

Z+ = z+(l + ^ 

V U r 

and linearizing under the assumption that u' TL <C u T gives 

ju*(z + )\ K L (t) 



' p ~ "MHM 



I + ( i + z+ d K ms (z + ) \ <{t) 

U*rms( Z+ ) dz + J U T 



+ [U*(z + ) + z> 



dz + J u T 
(12) 



From the relationship between r and u T after linearization one gets 

T U T 

From Section [2] we also have 

-^ = au' + OL . 



Hence, 



<^ ~ -au' + OL . (13) 



u T 2 

With this substitution, comparing (Q and f fT2|) gives 

a u (z + ) = ~(u*(z+) + z + ^^y (14) 

and 

^'i^ ^y )"' < i5) 
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Figure 1: Comparison of the theoretical expression of, (a) a u , and (b) f3 u , 
respectively from equations ( I14p and (j!5p with values found experimentally 
in Mathis et al. (2011). 



The universal mean profile U*(z + ) and the random mean square u* ms (z + ) 
of the universal fluctuation u*(. . . ,z + ) should, in principle, be adjusted so 
as to fit experimental results better over a range of Reynolds numbers. The 
most appropriate way would be of course to find the universal velocity u* 
by fitting (j5j) directly to experimental data, without introducing fluctuations 
and without linearization. However, due to the linearisation the results are 
approximate in any case. For the first check of the theory we took U* and u* ms 
to be equal to the mean profile in the calibration experiment and the random 
mean square of the universal signal of Mathis et al. (2011) respectively, and 
we used a = 0.0898 from the same source. Comparison of (fH|) and (TT5l) 
with the experimental results^] described in Mathis et al. (2011) is shown in 
Figures [TJ^a) and[T|6). 



4 Skin friction superposition coefficient 

So far, the skin friction superposition coefficient a was considered as a con- 
stant measured in experiment independently of all other quantities. Given 
its value, our theory established the relationship between the universal mean 
velocity profile and universal random-mean-square velocity fluctuation pro- 
file and the functions a u (z + ) and f3 u (z + ). For the purposes of comparisons 
in Figures Ufa) and QJ6) the experimentally measured value of a was used. 

One further step can be made if, in the spirit of the classical method of 
matched asymptotic expansions, one recognises that the logarithmic layer 
is, in fact, a part of the inner distinguished limito Because of that, the 
relationship f )12p . which is valid in the inner limit should also be valid in 
the logarithmic layer. In MHM the large scale signal u'q L was chosen to 
be the low-pass filtered velocity u't at a position in the logarithmic layer 
at z + = Zq L = 3.9Re7 2 . Applying the low-pass filter to ( TT2|) removes the 
fast-varying term. Substituting z + = Zq L and using ( TT3|) gives, 



■"■Note that the definition of a (a u ) in MHM should be corrected to include the ratio 
of random mean squares and should be a = max R(u~^,u'Q L )(u~^) rm .«/(u'oi)r.m.s' The 
plots were built with this correction. 

2 It is also the part of the outer distinguished limit and, as such, is the region of overlap 
of the two distinguished limits. 
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Figure 2: Dependence of a with the wall-normal location of Zq L . 
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U*(z + ) + z 
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Cancelling out u' OL one obtains 



a = 2/ I U*(z+) + z 



*(-y+' 



_dU*{z 
dz~ 



Since Zq L is inside the logarithmi c layer, one can take U*(z + ) = Mnz + + 
B, with k = 0.384 and B = 4.17 flMonkewitz. Chauhan. fc Nagib I IJOOrl ). 
This gives 

a = 2/ (-\nz+ L + B + -) . (16) 



K 



K 



Figure [2] shows the comparison for a calculated from ffTBI) along with 
values obtained from experiments and DNS data. It is seen that the value of 
a derived here depends on the choice of the wall-normal location Zq L at which 
the large-scales signal is taken. This result is consistent with MHM since the 
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coefficient a is proportional to the coefficient correlation max{R(u T ,uoL,)}, 
and is thus known to depend on the location of the outer-probe (Brown & 
Thomas; Marusic & Heuer 2007). 



5 Concluding remarks 

In this paper we have investigated the simplest theoretical description of the 
phenomenon of modulation of near-wall turbulence by large scale structures 
associated with the logarithmic layer. In this description the effect of large 
scale structures is limited to a quasi-steady response of the near- wall turbu- 
lence to slow large-scale fluctuations of the skin friction. The most natural 
and compact form of expressing this mechanism is given by (j2]) and (J5J), 
which are simply the usual Reynolds-number-independent representation of 
total skin friction and velocity in wall variables but with mean skin fric- 
tion replaced by large-scale low-pass-filtered fluctuating skin friction. This 
mechanism manifests itself both in amplitude modulation and superposition, 
introduced in previous works, and frequency and wavenumber modulation 
which has not been considered before. When rewritten in terms of fluctu- 
ations and linearised with an assumption that the large-scale fluctuation of 
the skin friction is small as compared to the mean skin friction, our formu- 
lae coincide with the formulae proposed by MHM. However, an additional 
feature here is that in the proposed theory the skin friction superposition 
and modulation constants a and (3 and the analogous functions a u (z + ) and 
fi u {z + ) for the velocity distributions are expressed explicitly via the univer- 
sal mean velocity profile, and the universal random mean square fluctuation 
velocity profile. In MHM these were obtained purely from experiments. The 
predicted values of a and (3 agree reasonably well with the values obtained 
experimentally. 

The agreement for a u (z + ) and f3 u (z + ) is also encouraging, especially if one 
takes into account the error introduced by the linearization. The linearization 
was needed, of course, only in order to be able to make the comparisons. It is 
more appropriate and might be easier to use the exact theoretical formulae. 
As far as the comparisons are concerned, one should also note that the theory 
is applicable only in the inner layer, including the log layer, and that the 
theoretical results expressed in wall units are independent of the Reynolds 
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number. The way (3 U was denned in experiments (MHM) was such that 
/3 u (z + = 3.9Rey 2 ) = 0, that is /3 U is dependent on the Reynolds number 
at least near the largest values of its argument, which, therefore, is outside 
the region of applicability of the theory. Overall, the comparisons strongly 
support the theory. 

It is of interest to note the simple physical mechanism that leads to the 
reduction, and even possible change of sign, of the amplitude modulation 
of the turbulence kinetic energy of small-scale fluctuations as the distance 
from the wall increases. An increase in the large-scale velocity has two ef- 
fects. First, it increases the skin friction, and that leads to an increase in the 
turbulence kinetic energy of small-scale fluctuations. Hence, the first effect 
creates a positive correlation between the large-scale velocity and small-scale 
turbulence intensity. Second, an increase in the skin friction leads to a de- 
crease in the thickness of the inner region. As a result, at a fixed physical 
distance from the wall the value of z + increases. If this point happens to be 
in the region where the turbulence kinetic energy of small-scale fluctuations 
decreases with z + , this second effect will decrease the degree of amplitude 
modulation, thus contributing negatively to the correlation. More accurate 
information on the universal distribution of the turbulence kinetic energy of 
small-scale fluctuations would be needed in order to confirm if this mech- 
anism can indeed lead to negative modulation at some distances from the 
wall. 

The modulation theory makes it possible to extrapolate the results of ob- 
servations at moderate Reynolds numbers to much higher values of Reynolds 
number. For example, it is well known that riblets reduce turbulent skin 
friction. The drag-reducing effect is strongest when riblet dimensions ex- 
pressed in wall units attain an optimal value. At high Reynolds number, 
however, the large-scale structures will result in slow variations of the skin 
friction, and, hence, this implies that the effective riblet dimension in wall 
units cannot remain optimal at all times, thus leading to the decrease of 
the effectiveness of riblets at higher Reynolds numbers. By predicting the 
magnitude of the large-scale struct ures, say, using the experimental fit of 



Alfredsson. Orlii. fc Segalini I (120 ll[ ). and then using the known dependence 



of the drag reduction effect on rible t geometry for moderate Reynolds number 



( iGarcia-Mayoral fc Jimenez 1 120111 ) and the formulae derived in the present 



study, it is possible to estimate quantitatively the change in the drag reduc- 
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tion due to an increase in Reynolds number. The same idea applies to drag 
reduction by wall oscillations. Of course, this implies that the modulation 
effect, confirmed by comparisons in the present work only for the flow past 
a flat solid wall, will also remain the only effect of Reynolds number in flows 
with drag reduction. These considerations show, however, the importance of 
further studies of the modulation effect. 
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